Electronic structure of relativistic Mott insulator L^RhOs 
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Motivated by studies of coexisting electron correlation and spin-orbit coupling effect in Na2lr03 
and a recent experiment of its 4d analogue Li2Rh03, we performed first-principles calculations of 
the rhodium oxide compound. The experimentally observed ground state of Li2Rh03 can be recov- 
ered only if both spin-orbit coupling and on-site Coulomb interaction are taken into consideration. 
Within the proper U range for 4d-orbitals (2 ^ U ^ 4 eV), the ground state of Li2Rh03 could 
be either zigzag-AFM or stripy-AFM, both yielding energy gap close to experimental observation. 
Furthermore, the total energy differences between the competing magnetic phases are ^ 3 meV/Rh 
within 2 ^ U ^ 4 eV, manifesting strong magnetic frustration in the compound. Finally, the 
phase energy of Li2Rh03 cannot be fitted with the two-dimensional Heisenberg-Kitaev model in- 
volving only the nearest neighbor interactions, and we propose that inter-layer interactions may be 
responsible for the discrepancy. 
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Electronic structure of transition metal oxides has been 
a hot topic in the condensed matter community for 
decades. On one hand, the 3d-transition metal oxides are 
considered as typical strongly correlated systems whose 
on-site electron-electron interactions U are large enough 
to forbid the hopping of electrons among lattice sites, 
leading to the Mott behavior. On the other hand, as the 
5d-orbitals are much more extended than the 3d-orbitals, 
the correlation effect for 5d-orbitals are small and negli- 
gible, while its spin-orbit coupling (SOC) effect becomes 
prominent. 

The above argument, however, is challenged by stud- 
ies of 5d- metal oxide Na2lr03 1 10 . As a 5d metal ox- 
ide, Na2lr03 is an ideal candidate for realizing the 
Heisenberg-Kitaev model 



#hk= [(i-^Si 



S j -2aS?S] 



(1) 



<*,J> 



where 7 = x,y,z denotes the three different types of links 
in the hexagonal lattices, and a indicates the interpola- 
tion between the Heisenberg term and the Kitaev term, 
with a = and a = 1 recovering the original Heisenberg 
and Kitaev model, respectively. It is believed that a is 
relatively large for 5d electron systems due to strong SOC 
effect. As pointed out by J. Reuther et al., the ground 
state of the Heisenberg-Kitaev model is Neel antiferro- 
magnetic (AFM) for < a < 0.4, and stripy AFM for 
0.4 $C a < 0.8. When the Kitaev term becomes domi- 
nant ly large where 0.8 ^ a < 1, the system will enter 
a spin-liquid phase 11 . However, it was later discovered 
in both first-principles calculations and experiments that 
the ground state for Na2lr03 is indeed the zigzag AFM 
phase, and the material is a relativistic Mott insulator—^. 
To resolve the apparent discrepancy, it was proposed that 
a next-nearest-neighbor Heisenberg coupling term should 
be added into the original Heisenberg-Kitaev model 7 . On 
the contrary, I. I. Mazin et al questioned the application 



of localized spin models in the system, and proposed a 
quasi- molecular-orbital approach^. 
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FIG. 1: (a) Crystal structure of Li 2 Rh03- The blue large 
atoms are Rh; red smaller atoms are Li; and yellow smallest 
atoms are O. The hexagonal lattice of Rh is also shown, (b) 
The first irreducible Brillioun zone and the high symmetry 
points appeared in this paper, (c-e) The Neel- AFM, zigzag- 
AFM, and stripy-AFM configurations for hexagonal lattices. 
The red and blue spheres indicates spin up and down Rh 
atoms, respectively. 

Recently, 4d-transition metal oxide Li2Rh03 was re- 
ported to be in close proximity to spin glassy phased. 
As the 4d systems usually have moderate Coulomb in- 
teraction and SOC in between 3d and 5d systems, this 
4d compound should be very interesting not only because 
of its structural analogue to Na2Ir03, but also due to the 
more delicate interplay between the Coulomb correlation 
and SOC. It was determined experimentally that no long- 
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-134.5490 


-134.6043 
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-134.6441 
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-135.2827 


-135.4302 
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-135.4610 
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-131.0182 
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-131.3295 
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TABLE I: Total energies (per unit cell, four Ru atoms) of different Li2Rh03 magnetic phases. Columns 2-6 are for the 
calculations without SOC; while columns 7-11 are for the calculations with SOC. The lowest phase energy at different U with 
and without SOC are indicated with italian and bold fonts, respectively. 
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FIG. 2: (a-b): Band structure of NM phase without (red solid lines) and with (blue dashed lines) SOC. (c-d): Total and partial 
DOS of NM phase with SOC. (a) and (c): U=0 eV. (b) and (d): [7=4 eV. 



range magnetic order was observed in Li2Rh03 down to 
0.5 K, and spin- freezing temperature was measured to 
be around 6 K. Fitting of the resistivity data suggested 
that the compound was a semiconductor with a narrow 
band gap of ~78 meV. It is worthy noting that anti-site 
disorder effect between Li + and Rh 4+ ions is inevitable 
in the experiment, therefore it is yet to be determined 
whether the spin-glassy feature is intrinsic or due to the 
anti-site defects. The experiment therefore raises several 
interesting questions: 1. Is the narrow band gap due to 
electron correlation or spin-orbit coupling effect? 2. How 
close is the defect-free Li2Rh03 to the spin-glassy phase? 
3. How relevant is the system to the Heisenberg-Kitaev 
model? 

In this article, we present our latest first principles 
study of the 4d-transition metal oxide Li2Rh03. We cal- 
culated and compared the total energy of its magnetic 



ordered phases with and without SOC effect as well as 
for different on-site interaction U . We analyzed its band 
structure and density of states (DOS) at different C/'s 
with and without SOC. We conclude that, both the SOC 
and on-site interaction U are important in the compound, 
and the energy gap is indeed a correlation effect. More- 
over, the magnetism in Li2Rh03 is extremely frustrated, 
and thus the system may be very close to spin disordered 
phase. Finally, attempts to fit the total energy data with 
Heisenberg-Kitaev model failed, suggesting that interac- 
tions beyond original Heisenberg-Kitaev model may be 
required to understand the magnetism of the system. 

To study the magnetic and electronic properties of 
Li2Rh03, we performed density functional based first 
principles calculations. In particular, we used the plane- 
wave basis and projected augmented wave method— as 
implemented in the Vienna Abinitio Simulation Package 
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(VASP)±ii£. The Perdew, Burke, and ErnzerhofT flavor 
(PBE)i£ of generalized gradient approximation (GGA) 
were chosen to be the exchange-correlation functional. 
To ensure the convergence of total energy to 1 meV/cell, 
a high energy cut-off of 540 eV to the plane wave basis 
was chosen; while a dense 8x5x8 T-centered K-grid 
was used to perform the Brillouin zone integration. The 
lattice constants and internal atomic positions were opti- 
mized so that the forces on individual atoms are smaller 
than 1 meV/A and internal stress less than 0.04 kbar. 
The density of states (DOS) were calculated using a 
16 x 10 x 16 T-centered K-grid and tetrahedron method. 

Solution of the Heisenberg-Kitaev model suggested 
three possible AFM ground states for hexagonal 
lattice o 1Q i 17 i 18 , which are also used to study the mag- 
netic phases of Na 2 Ir0 3 718 (FIG. |l(c)||l(e)| ). In our cur- 
rent study, we performed calculations for all three AFM 
states, as well as ferromagnetic (FM) state and non- 
magnetic (NM) state. The total energies of different 
magnetic phases are listed in TABLE U 

Our calculations demonstrate that the SOC effect do 
play an essential role in determining the ground state 
magnetism. The FM phase is always the ground state 
if SOC is not included in the calculation, although the 
magnetic frustration is greatly enhanced with respect to 
increased ?7, as indicated by the reduction of the total 
energy difference between the FM phase and the second 
lowest phase (from ~32 meV/Rh at U = eV to ~2 
meV/Rh at U = 4 eV). With the SOC being correctly 
accounted for, the ground state of Li2Rh03 starts from 
the FM phase at U = eV to the zigzag AFM phase 
at U = 2 eV, and eventually becomes stripy AFM phase 
at U = 4 eV. It is worthy noting that once the AFM 
ground state is established, the total energy difference 
between the ground state and second lowest phase does 
not exceed 12 meV/cell (or 3 meV/Rh) (11 meV/cell at 
U = 2 eV, 4 meV/cell at 17 = 3 eV, and 7 meV/cell at 
U = 4 eV). By comparison, the total energy difference 
between the NM phase and the highest magnetic phase is 
at least one order of magnitude larger (465 meV/cell at 
U = 2 eV, 850 meV/cell at U = 3 eV, 1253 meV/cell at 
U = 4 eV). Such small energy difference between mag- 
netic phases compared with the large magnetic energy 
manifests strong frustrations of magnetic interactions, 
and that the system may be very close to spin-glass or 
spin-liquid phases. 

We now turn to the electronic structure of Li2Rh03. 
Firstly, we compare the band structure of NM Li2Rh03 
calculated without (Fig. |2(a)[ ) and with (Fig. |2(b)[ ) SOC. 
Moderate band splitting due to SOC are normally ex- 
pected for 4d- transit ion metals. However, for the cur- 
rent system, the major SOC effect is the band reordering 
around T from Ep-0.5 eV to Ef+0.2 eV and from Ep- 
1.5 eV to E F -0.& eV. Nevertheless, a ~ 86 meV band 
splitting due to SOC can still be identified around Y at 
~ Ep — 1.0 eV, which signals the strength of the SOC. 



Li2Rh03 without and with SOC, respectively. Evidently, 
the electron states near Fermi level are contributed by the 
Rh-4d orbitals and 0-2p orbitals, which hybridize signif- 
icantly over a wide energy range. The t^g orbitals of 
rhodium due to the octahedral crystal field dominate the 
4d contributions. At the NM phase, the system remains 
metallic even LDA+Z7+SOC is employed with U up to 4 
eV. 
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Fig. |2(c) and 2(d) show the total density of states (DOS) 
as well as the projected density of states (PDOS) of NM 
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FIG. 3: Band structure of (a) zigzag and (b) stripy phases 
without considering SOC effect at different U level. The larger 
panel in each subfigure shows the band structure at U = eV; 
while the rest four smaller panels from left to right correspond 
to U =1, 2, 3, 4 eV, respectively. To save space, only the band 
structure from M' to Y is shown for U ^ cases. 



Secondly, we examine the electronic structure of 
Li2Rh03 at the zigzag- and stripy- AFM phases Al- 
though the AFM phases cannot be stabilized unless SOC 
is turned on, it is instructive to compare the electronic 
structure at different on-site U without SOC first (FIG. 
[3]). The large on-site energy pushes away the bands near 
the Fermi level, causing significant reduction of electron 
states at Ep. However, the electron-electron correlation 
alone cannot drive the system insulating, since the band 
gap is not opened even with U as high as 4 eV. With SOC, 
the system is also metallic at U = eV, although signif- 
icant band reordering effect similar to the one observed 
in the NM case is also present in both AFM cases (FIG. 
4(a)| and |4(b)[ ). However, a ~10 meV energy gap immedi- 
ately opens even with a small on-site energy U = 1 eV in 
the zigzag- AFM phase. The gap grows to ~75 meV and 
~130 meV, if the on-site energy U is increased to 2 eV 
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FIG. 4: Band structure of (a)(c) zigzag and (b)(d) stripy phases with SOC effect at different U level. 



and 3 eV, respectively. Noticing that the experimentally 
observed E g is 78 meV, which is very close to the cal- 
culated value for zigzag- AFM phase at U = 2 eV (FIG. 
4(c)] ). The magnetic moment is found to be — 0.4 /i#/Rh, 
twice the value for the one in the zigzag- AFM phase of 
Na2lr03 3 . For the stripy- AFM phase, a relatively larger 
U = 2 eV is required to open an diminishing energy gap 
of 3 meV. The gap expands to —62 meV and 117 meV 
for U = 3 eV and U = 4 eV, respectively. The mag- 
netic moment is —0.2 /i#/Rh in this phase, similar to the 
value for the one in the stripy- AFM phase of Na2lr03 3 . 
As the zigzag- AFM and stripy- AFM phases are energet- 
ically very close around U = 3 eV (— 1 meV/Rh), one 
cannot deduce the ground state magnetic ordering from 
the comparison between the calculation and experiment. 
Once again, the two phases severely competes with each 
other, and fluctuations (quantum or thermal) can easily 
drive the system into disordered spin states. Neverthe- 
less, in each case the system is non-metallic with a small 
gap of 60-120 meV. 

We want to make some further remarks regarding the 
obtained results. The calculated band structure in all 
magnetic phases appear to be quite three-dimensional. 
In fact, the inter-layer Rh-Rh distances are —5.1 A still 
comparable with the intra-layer Rh-Rh distance of 3.03 
A. Therefore, it is questionable whether two-dimensional 
Heisenberg-Kitaev model is applicable to such system. 
Indeed, when we attempted to fit the phase energy data 
to the classical results of the Heisenberg-Kitaev modelii, 
the fitting failed for all U, as the fitted J and a result 



in a wrong phase energy order. Therefore, it is possible 
that interactions beyond the original Heisenberg-Kitaev 
model is required to describe the system. Considering 
the three-dimensional feature of its electronic structure, 
we propose it is likely to be the inter-layer coupling J c . 

In conclusion, we have performed first-principles study 
of Li2Rh03, comparing the electronic structure and ener- 
getics of different magnetic phases with or without SOC 
effect at U ranging from to 4 eV. Both SOC and on-site 
electron correlation are crucial in determining the elec- 
tronic structure of the compound. Within the proper U 
range, the compound is semiconducting with a small gap 
of —60 to — 120 meV; and the energy difference between 
the competing magnetic phases are extremely small, sug- 
gesting highly frustrated magnetism in the system and 
that the system may be in close proximity to disordered 
spin state. Since the phase energies do not fit well with 
classical Heisenberg-Kitaev model results and its elec- 
tronic structure is highly three-dimensional, we propose 
that inter-layer coupling J c may be required to explain 
the behavior of the system. 
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